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Abstract. 

Context At least a fraction of the atomic hydrogen in spiral galaxies is suspected to be the result of molecular hydrogen which 
■^j- ' has been dissociated by radiation from massive stars. 

l/^ , Aims In this paper, we extend our earlier set of data from a small region of the Western spiral arm of M8 1 with CO observations 

' in order to study the interplay between the radiation field and the molecular and atomic hydrogen. 

Methods We report CO(1-0) observations with the Nobeyama 45 m dish and the Owens Valley interferometer array of selected 
regions in the Western spiral arm of M81. 

Results From our Nobeyama data, we detect CO(1-0) emission at several locations, coinciding spatially with H i features near a 
far-UV source. The levels and widths of the detected CO profiles are consistent with the CO(1-0) emission that can be expected 
from several large photo-dissociation regions with typical sizes of some 50 x 150pc located within our telescope beam. We do 
CIh' not detect emission at other pointings, even though several of those are near far-UV sources and accompanied by bright H i. 

JL , This non-detection is likely a consequence of the marginal area filling factor of photo-dissociation regions in our observations. 

We detect no emission in our Owens Valley data, consistent with the low intensity of the CO emission detected in that field by 
' the Nobeyama dish. 

Conclusions We explain the lack of CO(1-0) emission at positions farther from far-UV sources as a consequence of insufficient 

heating and excitation of the molecular gas at these positions, rather than as an absence of molecular hydrogen. 
K" ■ 

X 

5^ ' Key words, galaxies: individual: M81 - galaxies: ISM - ISM: molecules - ISM: clouds - radio lines: galaxies 
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1. Introduction 

Morphological evidence has been presented on several occa- 
sions over the past twenty years that a non-negligible fraction 
of the atomic gas found in the arms of spiral galaxies may be 
produced by the desttuction of molecular hydrogen in photo- 
dissociation regions (PDRs) surrounding young massive stars. 
This view, which was initially developed by comparing the rel- 
ative placement of dust lanes, H n regions, and H i density en- 
hancements in parts of M83 (Allen, Atherton & Tilanus 1985, 
1986) and M51 (Tilanus & Allen 1987), was extensively tested 
by Allen et al. (1997, hereafter AKBS97) in M81, based on the 
relative positions of a large number of far-UV (FUV), Ha, and 
Hi sources. AKBS97 interpreted the morphological informa- 
tion as implying that a significant part of the H i in the spiral 
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arms of M8 1 is a product of star formation rather than a pre- 
cursor to it. 

Nevertheless, by conventional measures, the amount of 
molecular hydrogen in M81 is not extensive. For example, CO 
observations from selected regions of M81, the locations of 
some of which are indicated in Fig. ^ have been reported by 
Brouillet, Baudry & Combes (1988), Sage & Westpfahl (1991) 
and Brouillet et al. (1991, hereafter B91), who find relatively 
weak CO(1-0) emission from the spiral arms on the basis of 
relatively low resolution single-dish observations. Taylor & 
Wilson (1998) reported the detection of resolved emission from 
three molecular clouds in one of those (Southern) regions with 
diameters of ~100pc, and Brouillet et al. (1998, hereafter B98) 
describe several clouds that they found interferometrically in 
one of the (Northern) regions where they had previously de- 
tected CO with a single dish. These data indicate that either the 
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Fig. 1. FUV contours overlaid on a grey-scale representa- 
tion of the Hi integrated intensity map of the Western 
part of the inner disk of M81. Overlaid are the fields of 
our NRO and OVRO observations, as well as the fields 
observed by Brouillet et al. (1998; labelled B98), and 
by Taylor & Wilson (1998; labelled TW98). Contour lev- 
els for the FUV are 2.3,3.1,4.6,6.2,7.7,9.2,13.9,18.5 and 
36.9 x 10~ 18 erg cirT 2 s _1 A -1 arcsec~ 2 . Grey levels for the Hi 
vary linearly from 1.5 (light) to 7.3 x 10 21 atoms cirT 2 (dark). 
FUV and H i data from AKBS97. 



earlier hypothesis regarding the origin of H i gas is incorrect for 
M81, or there must be a reservoir of (cold) molecular gas that 
has escaped detection. 

We have obtained new single-dish and interferometric mea- 
surements in the CO ( 1 -0) line from a region that extends below 
the southernmost field in the Western arm of M81 where B91 
detected CO emission, to determine whether evidence exists 
for a non-radiating molecular reservoir. The fields that we ob- 
served formed part of our previous study (AKBS97) and con- 
tain several candidate PDRs. Although CO gas spanning the H i 
velocity discontinuity caused by density wave-driven stream- 
ing of the gas (e.g., Visser 1980) was not detected from our 
interferometric observations, we obtained a number of single- 
dish detections near several Hn regions. We ascribe these de- 
tections to FUV excitation in the PDRs surrounding several 
young massive stars, but more spatially resolved data will be 
required at those sites to confirm this hypothesis. While our re- 
sults appear to corroborate the absence of molecular gas at or 
just beyond the spiral arm shock front, we argue that a popu- 



lation of large CO clouds similar to the cold Galactic PDR ob- 
served at l n * 216° by Williams & Maddalena (1996, hereafter 
WM96) in the plane of the outer Galaxy (i.e., G216-2.5) could 
exist at that location without being detected at the best synthe- 
sised spatial resolution or single-dish sensitivity achieved by 
our measurements. 

2. Empirical modelling of atomic and molecular 
volume densities 

Newly-formed massive stars produce copious amounts of FUV 
emission. This flux of FUV photons interacts with the surface 
of the parent molecular cloud and strongly affects its physics 
and chemistry by heating, dissociating, and exciting the gas 
(see Hollenbach & Tielens 1999 for an excellent review of the 
properties of such PDRs). When the presence of a PDR can be 
confirmed, e.g., from the expected "layered" morphology, use- 
ful information can be derived on the physical state of the gas 
from observations of the spectral lines emitted by the excited 
atoms and molecules in those regions. Of the many atomic con- 
stituents produced by photo-dissociation at and near the sur- 
faces of molecular clouds, atomic hydrogen (H i) is among the 
easiest to observe. 

The PDRs found previously in M81 by AKBS97 manifest 
themselves as bright "arcs" and "blisters" of enhanced 21 cm 
Hi emission, located near (and in some cases surrounding) 
sources of FUV continuum emission. Since the FUV emis- 
sion is not always accompanied by bright Ha emission, the 
sources of the dissociating photons must include (and may even 
be dominated by) B stars. This picture has been developed in 
more detail by Smith et al. (2000) who isolated 35 candidate 
PDRs in the giant spiral galaxy M101 and measured their Hi 
and FUV properties. A simple quantitative model for the Hi 
production in PDRs on the surfaces of molecular clouds was 
then used to determine the total volume density of the gas in 
the underlying molecular regions, which was found to lie in the 
range 30 - 1000 cirT 3 , without a radial dependence of the mean 
value. 

The production of Hi in PDRs by photodissociation has 
recently been examined in more detail by Allen, Heaton & 
Kaufman (2004, hereafter AHK04) who confirmed the appli- 
cability of a simple analytic formula (after some tuning of its 
parameters) initially proposed by Sternberg (1988). In addi- 
tion, AHK04 extended the picture to include CO(1-0) emis- 
sion, which is also produced in the same PDRs by the heat- 
ing of the ISM (mainly by grains via the photoelectric effect). 
AHK04 showed that observations of the H i column density and 
the CO(1-0) line intensity of the same PDR can be combined 
into a diagnostic which allows one to measure not only the den- 
sity of the underlying molecular cloud, but also the FUV flux 
incident upon it. They applied their approach to G216-2.5, and 
concluded that an incident FUV flux of ~0.8 (in units of the lo- 
cal FUV flux in the ISM near the sun) was bathing the surface 
of a large molecular cloud of density ~200cirT 3 , thereby cre- 
ating the observed levels of Hi column density and CO(1-0) 
surface brightness. 

Since the size of the Galactic PDR observed by WM96 
(100 - 200 pc) is of the same order as the linear resolution of 
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the earlier M81 observations analysed by AKBS97 (~150pc), 
we choose that object as a template for estimating the level 
of CO(l-O) emission we might expect in our targeted region 
of the spiral arms of M81. The total CO(l-O) emission from 
G216-2.5 is 7.1Kkms-' (as quoted by AHK04), and from 
Fig. 1 of WM96, we estimate that the total CO(l-O) pro- 
file width (FWHM) will be about lOkms -1 . From Fig. 2 in 
WM96, however, we estimate the extent of this CO(l-O) emis- 
sion to be ~50 x 150 pc 2 . Since our new single dish detec- 
tions in M81 make use of a millimetre radio telescope with 
an angular resolution of 16 arcsec » 280pc (1 arcsec = 17.5 pc 
at the 3.6 Mpc distance we adopted for M81; Freedman et al. 
1994), beam dilution would reduce the expected total emission 
to-O.SeKkms 1 . 

We have not attempted to compare this extrapolation with 
predictions from the simple model developed in AHK04 of 
the integrated CO( 1 -0) intensity expected near the shock front 
in M81 since that model includes photon heating only. The 
known correlation between the velocity-integrated CO inten- 
sity and the nonthermal radio continuum surface brightness 
in nearby galaxy disks (Adler, Allen, & Lo 1991) and in the 
Galaxy (Allen 1992) suggests that an enhanced level of cosmic 
rays also has something to do with heating the ISM, and there- 
fore with elevating the CO surface brightness. The AHK04 
model also utilizes the microturbulent velocity width of a sin- 
gle cloud rather than the macroturbulent width typifying an ex- 
tragalactic ensemble. Nevertheless, the faintness of nonthermal 
radio continuum emission from M8 1 indicates that it is a better 
choice for assessing the issues under investigation than other 
nonthermally-bright nearby galaxies such as M51, M83, and 
NGC 6946 because the reduced role of cosmic ray heating pro- 
vides a more controlled experiment. 



chopper- wheel method (Ulich & Haas 1976), yielding the an- 
tenna temperature T* K corrected for atmospheric attenuation. 
In this paper, we use the main beam brightness temperature 
r m b = T^/j] m \, as the intensity scale of the CO brightness tem- 
perature. The system noise temperature (SSB), including the at- 
mospheric effect and the antenna ohmic loss, was 500 - 800 K 
in 7*^, depending on the elevation of the source and the atmo- 
spheric conditions. We smoothed the spectra to a velocity res- 
olution of 10 kms" 1 . 

We observed 46 points in a region of about 70 x 90 arcsec. 
The separation of each point is 1 1 arcsec in right ascension and 
declination. The telescope pointing was checked and calibrated 
every few hours by observing the SiO maser emission of the 
late-type star R-UMa at 43 GHz. The absolute pointing accu- 
racy was better than 5 arcsec (peak value) throughout the ob- 
servations. 

At the data reduction stage, we checked all individual spec- 
tra, with integration times of 20 s, and eliminated those with 
bad baselines. The remaining spectra were coadded and fitted 
with a linear baseline. 

The region of interest is shown in detail in Fig.|2 where we 
show the FUV, Ha, and H i distributions at 9 arcsec resolution 
in the same format as used in AKBS97 for five larger regions 
of M81. Dust lanes (based on positions given by Kaufman, 
Elmegreen & Bash 1989a) are also indicated, in the middle 
panel. Figure [3] shows the NRO 45 m beam positions overlaid 
on the grey scale representation of the H i data, as shown in the 
right panel of Fig. [2] The beam positions designated with let- 
ters correspond to detections, as discussed in more detail below, 
and as tabulated in Tabled 



3. Observations, data handling, and results 

3.1. Nobeyama data 

Observations of the 12 CO(7 =1-0) emission were made dur- 
ing 1998 January 23 to 26 and March 27 to 31 with the 45 m 
telescope at the Nobeyama Radio Observatory (NRO). At the 
rest frequency of the 12 CO(1-0) transition (115.271204GHz), 
the half-power beam width (HPBW) is 16 arcsec. The aperture 
and main-beam efficiencies were T] a = 0.29 and rj m b = 0.40, 
respectively. 

We used the 2x2 Semiconductor - Insulator - 
Semiconductor (SIS) focal-plane array receiver which can si- 
multaneously observe four positions separated on the sky by 
34 arcsec each (Sunada et al. 1995). Differences in T* A be- 
tween the four beams due to the different sideband ratios 
were adjusted to be equal. We used the SIS single beam re- 
ceiver with a Martin-Puplett type SSB filter for an image side- 
band termination at the cryogenic temperature (4K) to ob- 
serve selected points effectively. As receiver backends we used 
2048 channel wide-band acousto-optical spectrometers (AOS). 
The frequency resolution and channel spacing were 250 kHz 
and 125 kHz, respectively. At 115 GHz the corresponding ve- 
locity resolution and velocity coverage are 0.65 km s~' and 
650 km s _1 . Calibration of the line intensity was made by the 



Pos. 


Offset 


/ T mb Sv 




dv 


*mb,peak 




(") 


(Kkms" 1 ) 


(kms-') 


(kms-') 


(mK) 


A 


5.5 -27.5 


2.32 


-6.3 


25.3 


90 ± 24 (3.7) 


B 


-5.5 -27.5 


2.73 


-4.0 


18.2 


145 + 23 (6.3) 


C 


5.5 -38.5 


2.52 


-14.6 


14.7 


187 + 24(7.8) 


D 


-5.5 -38.5 


1.86 


-6.6 


15.6 


155 + 32 (4.8) 


E 


-16.5 27.5 


1.20 


20.6 


7.7 


125 + 24 (5.2) 


F 


-38.5 -5.5 


1.40 


20.6 


14.3 


87 ± 25 (3.5) 


G 


-38.5 -27.5 


1.47 


10.5 


14.4 


85 ± 26 (3.2) 


H 


-5.5 -16.5 


1.07 


-9.1 


8.6 


75 ± 23 (3.3) 


I 


16.5 -49.5 


2.50 


-29.0 


19.5 


100 + 45 (2.2) 


J 


-5.5 -49.5 


3.19 


-38.8 


26.4 


90 + 46 (2.0) 



Table 1. CO(1-0) measurements obtained using the 45 m NRO 
dish, dv = 10 kms -1 , T m \, = 7^/0.40. The identification of the 
measurement, alphabetically as used elsewhere in this paper, 
and as a position offset from (0,0), which is RA= 09 h 50 m 44.0 s 
, Dec= 69° 17' 30" (B1950), is given in cols. 1 and 2; the 
integrated intensity in K kms -1 is in col. 3; the peak position 
and width of the profile, corresponding to velocity and velocity 
dispersion, are given in cols. 4 and 5; and the peak main beam 
temperature is in col. 6, where the error quoted is the rms of 
the spectrum (the number in parentheses is the ratio of the two, 
hence spectrum A is considered a 90/24 = 3.7<x detection). 
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Fig. 2. The region of the spiral arm in M8 1 studied in the present paper. Contours in all three panels are of the FUV emission, 

The contours are overlaid 
csecT 2 [dark]); 

Ha (middle panel; grey levels from 0.5 to 95 x 10~" ergcm^s -1 arcsec~ z ); and Hi (right panel, grey levels from 1.5 to 7.3 x 
10 21 atoms cm -2 ). The positions of our three most credible CO detections from our NRO observations (B, C, and D in Hig.|3j are 
indicated in each panel by open star symbols. Centre positions of the dust filaments as measured by Kaufman et al. (1989a) are 
indicated by the line segments drawn in the middle panel. 



and contour levels are 2.3,3.1,4.6,6.2,7.7,9.2, 13.9, 18.5 and 36.9 x 10" 18 ergcm^s" 1 A" 1 arcsec" 2 
on grey scale representations of the FUV (left panel, grey levels from 0.8 [light] to 15 x 10~ 18 erg cm~ 2 s _1 A 

v n ergcm _2 s _1 arcsec -2 " 1 



A sample of the final CO(1-0) spectra is shown in Fig. [4] At 
most of the 46 positions we did not detect any emission, which 
corresponds to upper limits ranging from 48 to 201 mK (3 <x). 
We did, however, detect CO emission at some of the positions, 
which are indicated by letters in Figs.[3]and[4] The Hi spectra 
shown in Fig.0 derived from a data set kindly made available 
by E. Brinks & F. Walter (2005, private communication), show 
that in all cases our detected CO emission is at the expected 
velocity. The most convincing of our detections are labelled B, 
C, and D (> 5<x detections in terms of T m ^ pea ^; see Table 0. 
While E also satisfies this 5cr detection criterion, it spans only 
a single channel on Fig. 4, and is therefore less convincing. 
The detections B, C, and D are clustered around a weak peak 
in the FUV emission at RA= 09 h 50 m 43 s , Dec= 69° 17' 10" 
(B 1950.0), and are located close to the Hi arc which almost 
completely surrounds this FUV peak. Around this cluster of de- 
tections there are four other, more tentative, detections, where 



an emission peak is seen at the expected velocity, but with a 
peak main beam brightness of only 2 - Act (A, H, I, J). At two 
further positions, labelled F and G, we detect emission with 
a peak main beam brightness of ~3cr. Peak antenna tempera- 
tures, r.m.s. noise levels, and results from Gaussian fits for all 
the pointings discussed above are summarized in Tabled 

3.2. Owens Valley data 

The millimetre array at the Owens Valley Radio Observatory 
(OVRO) was used on two separate occasions in its low spatial 
resolution mode to make CO(1-0) observations toward a re- 
gion in the Western arm of M81 (see Fig.^ which corresponds 
closely to the lxl arcmin box shown by B91 to the SE of 
their region Nl. The minimum and maximum baselines were 
15 m and 115 m during both sessions, corresponding approxi- 
mately to spatial scales of 620 and 80 pc, respectively, at the 
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H.A. (1950.0) 

Fig. 3. Positions of the NRO 45 m beam pointings, overlaid on 
the H i map as shown in the right panel of Fig.|2] Grey levels as 
in Fig. |2 Letters denote positions where CO was detected (see 
text); less secure detections are denoted by letters in parenthe- 
ses. 



adopted distance to M81. Tracks were obtained on 1998 May 
22-23 and 1998 November 5 comprising repetitive cycles of 
five separate 5 minute integrations on-target and two 5 minute 
integrations on the gain calibrator. Since the phase centre for 
the target observations was RA= 09 h 50 m 44.0 s , Dec= 69° 17' 
30" (B1950), our 62.9arcsec primary beam fell nearly com- 
pletely within the region observed with the NRO dish (Fig.^. 

The digital correlators at Owens Valley were configured 
to provide two slightly overlapping spectral bands in the up- 
per sideband, each with 32 channels of 1 MHz (2.6 kms 1 ) 
width. These bands were subsequently truncated and merged 
to provide 56 contiguous channels for each track. The data for 
each track were then merged to produce a dataset with 2480 
visibilities corresponding to 13.78 hrs of on-source integra- 
tion with a visibility-weighted average system temperature of 
>K. The resulting velocity coverage, centred at -30 kms -1 




-0.1 - 



-0.1 



-200 200 
^lsr (km/s) 



-200 200 
v lsr (km/s) 



Fig. 4. Final CO(1-0) spectra from NRO (drawn lines), 
smoothed to dv = 10 kms -1 velocity resolution. Letters in 
this figure correspond to positions of the telescope pointing as 
shown in Table 1 and Fig. [5] with the (0,0) offset corresponding 
to position RA= 09 h 50 m 44.0 s , Dec= 69° 17' 30" (B1950). 
They indicate spectra discussed in more detail in the paper; less 
secure detections are denoted by letters in parentheses. The dot- 
ted lines indicate H i spectra, on an arbitrary scale, at the same 
position. 



(relative to the local standard of rest), extended from +42.8 to 
-102.8 kms- 1 . 

All CO mapping was performed with the Astronomical 
Image Processing System (AIPS) package, using uniform 
weighting to preserve any extended faint emission in the field. 
The data were averaged in velocity to obtain a datacube with 
28 separate 5.2 kms -1 channels. The restoring beam was a 
Gaussian of width 3.7 x 3.6 arcsec (FWHM) at a position an- 
gle of 8.2 deg. The flux density noise in the dirty map was 
0.032 JybeairT 1 (uncorrected for primary beam attenuation) 
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which equates to a temperature noise of 221 mK in the main 
beam. In addition, the visibilities were tapered to also form a 
9x9 arcsec beam at a position angle of 46.6 deg for compari- 
son with the H i, Ha, and FUV datasets analysed previously by 
AKBS97. The noise in the latter (uniformly weighted, uncor- 
rected) dirty map was 0.055 Jy /beams 62.4 mK. Our resulting 
maps are 16.4 and 6.7 times more spatially resolved, respec- 
tively, than the single dish survey maps made in this region 
by B91. Since the latter had a sensitivity of « 20 mK and 
an on-target observation time of around 0.5 hr, equation C2 of 
AHK04 shows that our new 9 arcsec map has a flux density 
sensitivity that is approximately 14 times better than the previ- 
ous survey data. 

No clear evidence was obtained for CO emission to the 
limit of our sensitivity, neither in our 5.2kms~' channel maps 
prepared at 3.7 arcsec, nor in those at 9 arcsec resolution. 
While several features appear on each map, the number of 
features with positive flux density exceeding preset thresholds 
roughly equals the number of negative-valued threshold cross- 
ings. When the smoothed (9 arcsec) flux density maps were fil- 
tered at 3cr, no features repeated from one channel map to the 
next. 

To further quantify these results, we have derived an up- 
per limit to the CO(l-O) flux density from this region based 
on noise estimates made with the AIPS task imstat, which 
were derived by fitting Gaussian functions to pixel histograms. 
The corresponding 3<x upper limit to the flux density value 
is 163.8 mJybearrT 1 . Accordingly, the velocity-integrated 3<x 
upper limit in our field is (163.8 mJybeairT 1 ) x (2 x 10 6 Hz) 
= 3.276 x 10 5 (Jybeam _1 )Hz over the 2-channel span corre- 
sponding to our chosen velocity resolution. This corresponds 
to 0.85Jykms _1 or 1.03 K km s -1 for our 9 arcsec beam. 
Although several of our strongest detections at NRO had in- 
tegrated intensities greater than this level, only NRO pointings 
E, H, and part of B fell within our primary beam at Owens 
Valley. The integrated intensities from E and H are near our 
OVRO sensitivity limit, and B is on the periphery of the pri- 
mary OVRO beam. 



4. Discussion 

The results described in the previous section from the large dish 
at Nobeyama have extended the regions where CO has been 
found in the Western arm of M8 1 to locations that are farther 
down that feature than the southernmost field detected by B9 1 . 
While interesting in itself, the central issue in interpreting this 
result is: why did the emission stop there? The two most plausi- 
ble explanations are that the lack of CO(1-0) emission indicates 
the absence of measurable quantities of molecular gas in dark 
regions of the observed field, or that the field contains low den- 
sity gas that isn't heated sufficiently to reveal its presence radia- 
tively. The first interpretation is a consequence of the prevailing 
view that the CO(1-0) surface brightness and the column den- 
sity of molecular hydrogen are linearly related, independent 
of any effects caused by localised heating (Young & Scoville 
1991; Combes 1991). In this case, the H i ridge observed down- 
stream from the velocity shock front probably represents the 



compression of previously existing atomic gas (Kaufman et al. 
1989b). 

An alternative explanation that is consistent with our ob- 
servations is that the low flux of cosmic rays and the low FUV 
surface brightness are inadequate to sufficiently heat the ISM 
to shine in the CO(1-0) line in most of our field. A low cosmic 
ray flux is implied by the faint levels of nonthermal radio con- 
tinuum emission from M81, and low FUV levels are evident in 
many regions of Figure 2. Emission is observed in areas where 
heating by photons through the photoelectric effect, and by cos- 
mic rays through ionisation, leads to a rise in the kinetic tem- 
perature to a level sufficient, along with the density, to allow 
the radiating CO energy levels to be populated. 

To illustrate these points, we note that the cosmic ray inten- 
sity appropriate for the Galaxy in the neighbourhood of the Sun 
would heat giant molecular clouds (GMCs) that are just dense 
enough to have thermalized CO emission (n » 10 4 cirT 3 ) only 
to ~5 K (Tielens 2005, p. 352). This is much lower than is ob- 
served in localized regions of GMCs in the Galaxy. At densities 
below ~10 3 the CO level populations would be sub thermally 
excited, leading to even lower radiation temperatures. As justi- 
fication for comparing M8 1 with the Milky Way, we note that 
the radio continuum surface brightness, which can be used as 
a proxy for the cosmic ray flux, is similar in the disks of our 
Galaxy and M81, at a level that is an order of magnitude lower 
than that in M51 (Adler et al. 1991; Allen 1992). 

We can readily draw two conclusions from this discussion. 
First, regions showing high-temperature CO emission must be 
heated by a supplementary source in addition to cosmic rays 
- this source is most likely photoelectric heating of the gas 
caused by dust absorption of FUV photons, although in some 
cases H2 formation heating may play a role (Rollig et al. 2006). 
In Galactic GMCs, such sources can usually be easily identi- 
fied as nearby recently-formed massive O/B stars. Second, in 
the absence of such photons and a sufficient level of cosmic ray 
flux, the cooling in GMCs is so efficient that the kinetic temper- 
ature will quickly fall below 5 K, making the detection in any 
emission line tracer virtually impossible. This, we postulate, is 
what happens in M8 1 . 

The lack of CO emission at locations between the H 1 veloc- 
ity shock and the atomic hydrogen ridge, and its weak detection 
near several H 11 regions downstream from the shock, is consis- 
tent with the interpretation by Allen et al. (1985, 1986) and 
Tilanus & Allen (1987) of morphological data from M83 and 
M51, i.e., that the interstellar gas remains primarily molecular 
after passing through the shock until it dissociates in the neigh- 
bourhood of giant H n region complexes downstream. The spa- 
tial extent of our observations is more than adequate to search 
for such patterns. Since one arcsec corresponds to ~17.5pc at 
the 3.6 Mpc distance assumed for M81, our 63 arcsec OVRO 
field alone spans nearly 1.1 kpc. This is more than twice the 
separation claimed by Kaufman et al. (1989b) for the distance 
between the H 1 velocity shock and the atomic hydrogen ridge 
in the Western arm of the latter galaxy. 

The cold Galactic cloud described by WM96, G216-2.5, 
provides a useful template for what we might observe from 
low density gas clouds in M8 1 . Although AHK04 found that 
the most likely density for this cloud was 200cirT 3 , the un- 
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certainty in that determination was large enough to render 
the role of collisional excitation uncertain. Since the beam- 
diluted integrated intensity of this PDR at the distance of 
M81 would be ~0.9Kkms _1 when observed by the 16arcsec 
NRO beam, similarly excited clouds with similar profile widths 
(~10kms _1 , cf. Sect. 2) should be apparent in our observa- 
tions. The three firm detections found earlier (i.e., B, C, and 
D) are all somewhat brighter (by a factor of 2 - 3) and broader 
(by a factor of ~1.5) than this estimate, perhaps indicating the 
presence of multiple PDR structures at these positions. 

The absence of CO emission in our single beam fields that 
cover the region between the velocity front and the H i ridge 
sets an upper limit of ~0.8Kkms _1 for the integrated inten- 
sity of the gas from that region. If low density molecular gas 
does reside there, it is not excited sufficiently by the low FUV 
flux to be detected at the best sensitivity we achieved. The lack 
of emission at more strongly irradiated locations along the H i 
ridges in Fig. 3 could easily be caused by enhanced beam di- 
lution from structures that are smaller in linear size than the 
WM96 cloud. The WM96 cloud would have a beam-diluted 
integrated intensity ~2.7Kkms~ 1 in our 9arcsec synthesized 
OVRO beam at the distance of M8 1 . Although this exceeds the 
3cr sensitivity value determined for that array by a factor of 
~2.6, the pointings detected by the Nobeyama telescope only 
partially fell within the OVRO field. 

To test the consistency of our findings with earlier results, 
and the plausibility of using the WM96 PDR as a template for 
modelling putative undetected gas clouds in our region of the 
Western arm, we reinterpret the results obtained by B98 with 
the IRAM Plateau de Bure interferometer. The location of their 
field is shown in our Fig. 1, while the resulting CO(l-O) map 
appears in Fig. 1 of their paper. From the extent and maximum 
intensity of the structures shown in their figure, we estimate 
that the NRO 45 m beam would have reported a brightness of 
only ~0.4Kkms _1 for their region C3. The combination of 
regions CI and C2, which could fall into one of our beams, 
would have an integrated intensity of about 1 Kkms -1 . Thus, 
we would not have detected clouds in our area with the charac- 
teristics found by B98 for C3, while the combination of clouds 
CI and C2 would have yielded only a marginal detection. None 
of these clouds would have been detected at the 3<x level in our 
9 arcsec OVRO maps. 

When we compare the extrapolated properties of the 
G216-2.5 PDR with the earlier CO results from B98, it is clear 
that they have similar CO brightnesses (~0.8Kkms _1 for the 
latter) and sizes (50x1 00 pc for the WM96 cloud, and at most 
150 - 250pc for the B98 clouds). The general lack of CO de- 
tections in the area observed with the NRO dish is entirely con- 
sistent with this picture: we detect a few stronger peaks, but 
fail to detect CO at the majority of our pointings. This is ex- 
actly what we would expect to see if we had observed clouds 
in our region that were similar to those found by B98 in their 
field: a marginal detection from C 1 +C2, and no detections else- 
where in the field. We thus ascribe our lack of CO detections 
to being predominantly a filling factor effect — individual CO- 
emitting regions may be relatively bright, but in our 16 arcsec 
NRO beam they would be diluted to below the detection thresh- 
old. 



5. Conclusions 

We have detected CO(l-O) emission at several locations show- 
ing bright H i features near FUV sources in M8 1 . The levels 
and widths of the detected CO profiles are consistent with the 
CO(l-O) emission to be expected from several large photo- 
dissociation regions with typical sizes of ~50 x 150pc located 
within our NRO beam of size 280 pc. Non-detections, either 
at other pointings near FUV sources that contain bright H i, or 
in our synthesized maps, are likely to be a consequence of the 
marginal area filling factors of PDRs in our observations, and 
future surveys should be prepared to substantially increase the 
sensitivity and spatial resolution over what we have achieved 
here. The lack of CO(1-0) emission at positions farther from 
the sources of FUV is consistent with insufficient heating of the 
ISM, resulting in insufficient excitation of the CO molecules. 
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